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We apply the ab initio no-core shell model/resonating group method approach to calculate the 
cross sections of the ^H(d,n)*He and ^He(ci,p)*He fusion reactions. These are important reactions 
for the Big Bang nucleosynthesis and the future of energy generation on Earth. Starting from 
a selected similarity-transformed chiral nucleon-nucleon interaction that accurately describes two- 
nucleon data, we performed many-body calculations that predict the S-factor of both reactions. 
Virtual three-body breakup effects are obtained by including excited pseudostates of the deuteron 
in the calculation. Our results are in satisfactory agreement with experimental data and pave the 
way for microscopic investigations of polarization and electron screening effects, of the ='H(d,7)^He 
radiative capture and other reactions relevant to fusion research. 



The ^H((i,n)^He and '^He((i,p)'*He reactions are leading 
processes in the primordial formation of the very light 
elements (mass number, A < 7), affecting the predic- 
tions of Big Bang nuleosynthesis for light nucleus abun- 
dances With its low activation energy and high 
yield, ^H((i,n)^He is also the easiest reaction to achieve 
on Earth, and is pursued by research facilities directed 
toward developing fusion power by either magnetic {e.g. 
ITER [2) or inertial {e.g. NIF Q) confinement. The 
cross section for the d-f-'^H fusion is well known exper- 
imentally, while more uncertain is the situation for the 
branch of this reaction, '^H((i, 7)'^He, that is being con- 
sidered as a possible plasma diagnostics in modern fu- 
sion experiments @ . Larger uncertainties dominate also 
the ^IIe((i,p)^He reaction that is known for presenting 
considerable electron-screening effects at energies acces- 
sible by beam-target experiments. Here, the electrons 
bound to the target (usually a neutral atom or molecule) 
lead to enhanced values (increasingly with decreasing en- 
ergy) for the reaction-rate, effectively preventing direct 
access to the astrophysically relevant bare-nucleus cross 
section. Consensus on the physics mechanism behind 
this enhancement is not been reached yet ^ , largely be- 
cause of the difficulty of determining the absolute value of 
the bare cross section. Past theoretical investigations of 
these fusion reactions include various i?-matrix analyses 
of experimental data at higher energies as well as 
microscopic calculations with phenomenological interac- 
tions Ul, ll2|. However, in view of remaining experimen- 
tal challenges (some of which described above) and the 
large role played by theory in extracting the astrophysi- 
cally important information, it would be highly desirable 
to achieve a microscopic description of the '^H(d,n)^He 
and '^He((i,p)'*Hc fusion reactions that encompasses the 
dynamic of all five nucleons and is based on the fun- 
damental underlying physics: the realistic interactions 
among nucleons and the structure of the fusing nuclei. 

In this Letter, we present the first ab initio many-body 
calculation of the ^H(d,n)'*He and ^He((i,p)^He fusion re- 
actions starting from a nucleon-nucleon {NN) interaction 



that describes two-nuclcon properties with high accuracy. 
The present calculations are performed in the framework 
of the ab initio no-core shell model/resonating-group 
method (NCSM/RGM) [iMi, a unified approach to 
bound and scattering states of light nuclei. We use, in 
particular, the orthonormalizcd many-body wave func- 
tions {v being the channel index) 



(1) 



with inter-cluster antisymmctrizcr for the {A—a,a) par- 
tition Au, ccnter-of-mass separation fA-a,a^ and binary- 
cluster channel states 
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The inter-cluster relative-motion wave functions xi^'^ ('') 
satisfy the integral-differential coupled-channel equations 



Y,Jdr'r'^[f^-^-nM~^U,{r, 



j^Xu'{r') 



with bound- or scattering-state boundary conditions. 
Here H^^F {r,r') and M;^^,^ {r,r'), commonly referred to 
as integration kernels, are respectively the Hamiltonian 
and overlap (or norm) matrix elements over the antisym- 
metrized basis of Eq. ^ . They contain all nuclear struc- 
ture and antisymmetrization properties of the problem. 

In the present application we investigate reactions 
involving A—^ nucleons, characterized by a dcuteron- 
nucleus entrance and nucleon-nucleus exit channels [a=2 
and a=l in Eq. respectively]. The NCSM/RGM for- 
malism for an a=l projectile was presented in Ref. [isj 
(where the interest reader can find further details on the 
approach), while the deuteron projectile formalism was 
introduced in Ref. [l^, where we investigated the c?-^He 
system. To calculate the '^H(d,n)''He and '^He(c?,p)^He 
reactions, we had to address the additional contributions 
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TABLE I. Calculated g.s. energies and point-proton rms radii 
of ^H, ^H, ^'He, and ''He obtained by using the SRG-N=*LO 
NN potential with A=1.5 fm~^ are compared to the corre- 
sponding experimental values. The NCSM calculations were 
performed in an HO space with Ninax~12 and M}=14 MeV. 
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of matrix elements (two for the norm and 5 for the Hamil- 
tonian kernel, respectively) between the two mass parti- 
tions: (A— 1,1) and (A— 2,2). These extensions of the 
formalism will be described in detail in a separate paper. 

The input into Eq. ^ are: («) the nuclear Hamil- 
tonian, particularly the chiral N^LO NN potential of 
Ref. [16|, which we soften by a similarity renormalization 
group (SRG) transformation [l3,[l3 characterized by an 
evolution parameter A; and (ii) the eigenstates of the 
interacting nuclei, i.e. ^H, ^H, "^He and ^Hc, calculated 
within the NCSM In this first attempt of provid- 

ing an ab initio description of the d+'^H (d+'^He) fusion, 
we omit both the chiral and SRG-induced three-nucleon 
(NNN) force components of the Hamiltonian, and se- 
lect a value of the SRG parameter (A = 1.5 fm^^) for 
which we reproduce the experimental Q value of both 
reactions within 1%. While a complete (A- independent) 
calculation should include these terms (and efforts in this 
direction are under way) , we argue that this is a fair ap- 
proximation for the time being. Indeed, for these very 
light nuclei the initial attractive chiral NNN force can- 
cels to some extent with that induced by the SRG evo- 
lution of the potential, which is repulsive in this A 
range [20| . Further, even though the fusion proceeds at 
very low energies, the deformation and virtual breakup 
of the reacting nuclei cannot be disregarded, particularly 
for the weakly-bound deuteron. A proper treatment of 
deuteron-breakup effects requires the inclusion of three- 
body continuum states (neutron-proton-nueleus) and is 
very challenging. In this first fusion application we limit 
ourselves to binary-cluster channels and approximate vir- 
tual three-body breakup effects by discretizing the con- 
tinuum with excited deuteron pseudostates, strategy that 



proved successful in our d-'^He calculations 15] 

We start by discussing our results for the ground states 
(g.s.) of d, "^H, ^He and ^He, the energies and radii of 
which are compared to experiment in Table HI The soft 
NN interaction (SRG-N^LO with A=1.5 fm"!) and har- 
monic oscillator (HO) frequency (/if2=14 MeV ) ad opted 
are the same as in the d-^He study of Ref. [15|. En- 
ergy convergence (at < 20 keV level) of the present re- 



120 

90 
60 

M 30 
u 

T3 

w 

-30 
-60 
-90 



60 



M 30 



3/2 



1/2 



4 4 

n+ He n+ He 



-- rf+'H 

■ ■ ■ no coupling 




-30 





\^ — .^d+'He ^ d+' He 


■' ■/ 4_ 

■ 1 ■' ^ 

■ 

' A 


3/2 " ~ r??**^^'^;^?^-; 






■1' /' 

/ J' 


— n(p)+*}ie 
■ ■ ■ no coupling 




's " 

1/2 

(b) 







0.5 



1.5 



£,,„[MeV] 



FIG. 1. (Color online) Calculated elastic n-*He (a), d-^H 
and d-^He (b) phase shifts. The dashed (dotted) lines are ob- 
tained with (without) coupling of the n(p)-*He and ci-^H(^He) 
channels and all nuclei in their g.s. The full lines represent 
calculations that further couple channels with one ^Si— ^Di 
deuteron pseudostate. The SRG-N^LO NN potential with 
A=1.5 fm^^ and the HO space with Afmax=12 (Aniax=13 for 
the negative parity) and HQ— 14 MeV were used. 



suits is reached for an HO basis size of A',„ax=12, where 
we also find a weak frequency dependence in the range 
ll</if7<18 MeV. 

Next, we consider the elastic phase shifts for both en- 
trance and exit channels. In the past, we had already 
studied n(p)-''He scattering within the NCSM/RGM [11, 
14{ . Here, we extend those calculations by including the 
coupling to the d-^H (d-^Hc) channels. The impact of 



this coupling can be judged (in the n+ He case) from 
Fig. (Ha). Besides a slight shift of the P-wave resonances 
to lower energies, the most striking feature is the ap- 
pearance of a resonance in the "^0^/2 partial wave, just 
above the d-^H (d-'^He) threshold. The further inclu- 
sion of distortions of the deuteron via an ^H 
pseudostate (d*), enhances this resonance, while leaving 
the other partial waves mostly unaffected. In the d-^H 
(d-^He) elastic scattering phase shifts of Fig. [TJ6) we ob- 
serve a moderate impact of the coupling to the n-^He 
(p-^'He) channels on the S'-waves. As for the nucleon- 
''He case, the ^6*1/2 phase shifts are repulsive due to the 
Pauli blocking (which is treated exactly in our formal- 
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FIG. 2. (Color online) Calculated S-factors of the ^H(d,n)''He 
(a) and ■^He(d,p)''He (b) reactions compared to experimental 
data. Convergence with the number of deuteron pseudostates 
in '^Si—'^Di (d*) and "^02 {d'*) channels is shown. See also 
caption of Fig.[l]for details on interaction and HO space used. 



ism). There is, however, a resonance close to threshold in 
the ^^83/2 channel (where projectile and target spins are 
aligned) that is enhanced by distortions of the deuteron. 

Finally, from the scattering phaseshifts we obtain the 
■^H((i.n)'*Hc and ■^Hc((i,p)^He cross sections. The corre- 
sponding S-factors are compared to various data sets in 
panels (a) and (b) of Fig. ^ respectively. The deuteron 
deformation and its virtual breakup play a crucial role. 
We show in particular the dependence on the number 
of pseudostates in the ^Si-^Di (d*) and {d'*) 
channels, included in the calculation. Energies of these 
pseudostates can be found in Table II of Ref. [l^ ■ The 
S-factors increase dramatically with the number of pseu- 
dostates until convergence is reached for 9d* + 5d'* . Our 
calculation depends also on the size of the HO basis used 
to expand the eigenstates of the reacting nuclei as well 
as the localized parts of the integration kernels (see Eqs. 
(5), (6) and Sec. II. B. of Ref. [l^). As for the bound 
states, we find a satisfactory convergence [see Fig. |3](a)]. 
Before demonstrating this point in more detail, here we 
would like to discuss the comparison with data. 

The experimental position of the '^He((i,p)'^IIc S-factor 
maximum is well reproduced (within few tens of keV) 



in our calculations [Fig. [2] (&)]. Overall, the agreement 
with experiment is quite reasonable, except at very low 
energies where the beam-target data are enhanced by 
the electron screening. For the '^H((i,n)'*He S-factor, the 
absolute difference between theoretical and experimental 
peak positions (~ 10 keV) is of the same order of magni- 
tude found in the d-l-'^He case, however the relative dif- 
ference is much larger for such a low-energy resonance. 
As a consequence, the ■^II((i,n)^He S-factor maximum is 
somewhat underestimated in our calculations and, hence, 
the calculated S-factor underestimates the data below 
^ 70 keV. The inclusion of the NNN force (chiral and 
SRG-induced) into the calculation should provide closer 
agreement with experiment, although possibly it would 
require an NNN interaction accuracy beyond what is 
currently used in theoretical nuclear physics. However, 
an even more important consideration is the following. In 
obtaining the eigenstates of the reacting nuclei, we take 
into account Coulomb and isospin breaking of the NN 
interaction. At the same time, we perform isospin pro- 
jections when evaluating the NCSM/RGM kernels. It is 
therefore understandable that the splitting between the 
two peaks may become slightly underestimated in our 
calculations, so that it is hard to reproduce them equally 
well simultaneously and a certain amount of tuning of 
the nuclear interaction may be unavoidable. 

To reproduce the position of the '^H(c?,n)^He S-factor 
maximum, we performed additional calculations using 
SRG-N^LO NN potentials with a lower A. In Fig. [3l we 
show that using A=1.45 fm"-'^, we are able to reproduce 
the experimental position of the maximum (we find also 
a 0.6% variation of the calculated Q value, towards even 
closer agreement with the measured one). The theoret- 
ical S-factor is then in an overall better agreement with 
data, although it is slightly narrower and its peak is some- 
what overestimated [Fig. [3] (a)]. This calculation would 
suggest that some electron screening enhancement could 
be also present in the '^H(c?,n)^He measured S-factor be- 
low - 10 keV. 

Finally, the convergence of the calculation with respect 
to HO basis size and number of deuteron pseudostates is 
very similar for the two A values considered. In Fig.|3](a), 
we present the '^H(c?,n)*He S-factor dependence on the 
size of the HO basis for iVi„ax=8— 12. We find a satis- 
factory convergence and expect that an iVinax=14 calcu- 
lation, which is currently out of reach due to computa- 
tional reasons, would not be significantly different from 
the present results. Also, in Fig. [3] {b) we show the con- 
vergence of the and ^D^li phase shifts with the 
number of deuteron pseudostates in the vicinity of the 
3/2+ •^H(d,n)''He resonance. This picture is also interest- 
ing as it highlights how the '^H(c?,n)*He and '^He(c?,p)^He 
fusion processes proceed through the 5*3/2 resonance in 
the entrance channel and the ^I?3/2 resonance in the exit 
channel. The tensor interaction, which is automatically 
included in the accurate NN potentials we are using, is 
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FIG. 3. (Color online) Calculated S-factor of the ^H(d,n)'*He 
reaction compared to experimental data (a) and diagonal 
^^3/2 n-'^He and '^83/2 d-'^H phase shifts (6). Convergence 
with A'^max and the number of deuteron pseudostates in 
^Si-^Di (d*) and {d") channels are also shown in (a) 
and (b), respectively. The A'inax=8, 10, and 12 results con- 
tain 9d* plus 3,4, and 5d'*, respectively. The n-*He kinetic 
energy is shifted by the d— '^H threshold energy. The SRG- 
N^LO iViV potential with A=1.45 fm"^ and the HO frequency 
?lO=14 MeV were used. 



TABLE II. Calculated S-factors at zero energy compared to 
the 7?-matrix data evaluation of Ref. O. The NCSM/RGM 
calculations as described in Figs. [3] and [2] for ''lI(ci,n)^He and 
^He(d,p)'*He, respectively. 



indispensable for the reaction to take place. Unlike the 
5*3/2, the ^£'3/2 phase shift does not cross 90 degrees, 
remaining positive near the resonance. We note the sim- 
ilarity of our calculated phase shifts with those extracted 
from the data by using the single-level i?-matrix fit of 
Ref. [1]. In Table HH we summarize our S{0) values and 
compare them to the i?-matrix analysis of Ref. [l3l . 

In conclusion, we performed ab initio many-body cal- 
culations of the ■^H(c?,n)^He and ■^He((i,p)'*He fusion re- 
actions. Our results are promising and pave the way 
for microscopic investigations of polarization and elec- 
tron screening effects, of the ■^H((i,7)^He radiative cap- 
ture and other reactions relevant to the fusion research 
that are less well understood or hard to measure. Our cal- 
culations can be further improved by including additional 



five-body correlations, e.g., virtual breakup of ("^He). 
This can be best done by coupling the NCSM/RGM 
binary-cluster basis with the NCSM calculations for ^He 
(^Li) as outlined in Ref. [2l|. Virtual excitations of the 
deuteron should be treated by considering explicitly n-p- 
'^H('^He) three-cluster channels. The inclusion of NNN 
interactions, both chiral and SRG-induced [2^, is also 
desirable. Efforts in these directions are under way. 
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